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SECTION  I 
INTRODUCTION 

PURPOSE  AND  SCOPE 

Because  of  the  elastic  properties  of  most  acoustic  targets, 
the  temporal  echoes  and  frequency  characteristics  of  even  the  simplest 
targets  are  distinctive  but  confusingly  complex.  These  same  effects  can 
influence  acoustic  images  as  well.  Any  useful  target  classification 
algorithm,  via  temporal  or  image  processing,  should  be  based  upon  mean¬ 
ingful  target  scattering  characteristics.  Hence,  the  frequency  response, 
or  differential  scattering  cross-section,  of  targets  of  interest  need  to 
be  measured  and  understood  so  that  it  may  guide  pattern  recognition 
procedures  to  those  frequency  domains  which  present  the  most  consistent 
and  distinctive  features.  Once  these  features  have  been  identified,  a 
waveform  may  be  designed  to  maximize  the  desired  response. 

The  purpose  of  this  report  is  to  present  both  mathematical  and 
heuristic  descriptions  of  acoustic  scattering  which  will  be  useful  in 
the  interpretation  of  experimental  scattering  data  from  minelike  objects 
by  providing  physical  insight  toward  understanding  their  origin  and 
significance.  Subsequent  reports  will  deal  with  the  results  of  experi¬ 
mental  measurements  and  their  interpretation. 

This  report  is  organized  into  nine  sections.  In  addition  to  these 
introductory  remarks,  Section  I  contains  a  brief  description  of  previous 
work  in  acoustic  scattering  and  the  phenomenology  involved.  Section  II 
is  a  derivation  of  the  general  scattering  theory  from  the  wave  equation 
to  the  general  Kirchhoff  diffraction  equation.  Thi.';  derivation  is 
included  for  continuity  <.nd  completeness.  Section  III  is  concerned  with 
the  development  of  an  original  empirical  scattering  model  based  upon  a 
basis  set  of  plane  waves  and  target  surface  distribution  coefficients. 
Section  IV  deals  with  a  method  of  presenting  and  interpreting  scattering 
data  in  terms  of  real-space  and  wave  (or  Fourier)-space  functions. 

Section  V  presents  a  derivation  of  a  first  level  approximation  to  scat¬ 
tering  from  a  rigid,  finite  cylinder  from  arbitrary  aspect  angle,  using 
the  empirical  model  described  in  Section  III.  Section  VI  contains 
graphical  displays  of  the  scattering  function  derived  in  Section  V  as 
well  as  comments  on  the  characteristics  of  the  plotted  function.  Sec¬ 
tion  VII  considers  impedance  boundary  conditions  for  elastic  targets. 
Reflection  coefficients  for  the  finite  cylinder  are  derived  in  terms  of 
the  wave-space  and  real-space  parameters  of  the  empirical  model.  Sec¬ 
tion  VIII  contains  plots  of  the  finite,  impedance  cylinder  (aluminum  and 
brass)  in  water  as  well  as  a  comparison  of  that  for  an  aluminum  cylinder 
in  air  with  that  of  the  rigid  cylinder  (Section  VI).  Section  IX  concludes 
with  a  summary  of  the  empirical  model  and  its  utility. 
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PREVIOUS  WORK  AND  PHENOMENOLOGY 

Acoustic  scattering  measurements  for  simple  targets  have  been 
published  in  the  form  of  time  series  (temporal)  signals  for  short  pulses^ 
or  as  power  spectra  (frequency  domain)  of  the  scattered  signal. 

The  backscattering  experimental  configuration  has  been  the  principal 
method.  Spheres  and  "infinite"  cylinders,  from  beam  aspect,  have  been 
the  principal  targets  because  their  symmetry  permits  a  minimum  number  of 
measurements  because  of  the  absence  of  aspect  dependence.  Some  early 
work^’®  shows  target  directivity  patterns  in  the  form  of  polar  plots  at 
discrete  frequencies.  In  these  investigations,  the  acoustic  projector 
and  target  were  held  fixed  while  the  receiving  hydrophone  was  rotated 
about  the  axis  of  symmetry  of  the  target.  Hence  these  results  show  a 
variation  of  target  response  due  to  changing  the  experimental  geometry 
(bistatic  angle)  but  not  because  of  the  target's  geometry  or  aspect 
dependence.  Little  work  has  been  reported  which  combines  the  variation 
of  the  target  response  with  respect  to  frequency  and  target  aspect  for 
either  backscattering  or  bistatic  angle  experimental  configurations. 

Dunsiger^  has  examined  the  aspect  dependence  of  a  variety  of  simple 
shapes  using  short  pulses  at  high  frequency,,  where  frequency  dependence 
of  the  target  response  was  assumed  to  be  unimportant.  The  resulting 


^Barnard,  G.  R.  and  McKinney,  C.  M.,  "Scattering  of  Acoustic  Energy  by 
Solid  and  Air  Filled  Cylinders  in  Water,"  J.  Acoustic.  Soc.  Am.  33, 

226-238  (1961). 

^Hickling,  R.,  "Analysis  of  Echoes  from  a  Solid  Elastic  Sphere  in  Water," 

J.  Acoust.  Soc.  Am.  34,  1582-1592  (1962). 

^Diercks,  K.  J.  and  Hickling,  R.,  "Echoes  from  Hollow  Aluminum  Spheres  in 
Water,"  J.  Acoust.  Soc.  Am.  4l,  380-393  (1966). 

^Shirley,  D.  J.  and  Diercks,  K.  J.,  "Analysis  of  the  Frequency  Response  of 
Simple  Geometric  Targets,"  ARL-TM-69-21,  The  University  of  Texas  at  Austin, 
45  pp.  (1970). 

®Neubauer,  W.  G.,  Vogt,  R.  H. ,  and  Dragonette,  L.  R.,  "Acoustic  Reflection 
from  Elastic  Spheres.  I.  Steady-State  Signals,"  J.  Acoust.  Soc.  Am.  55, 
1123-1129  (1974). 

®Faran,  J.  J.,  "Sound  Scattering  by  Solid  Cylinders  and  Spheres,"  J.  Acoust. 
Soc.  Am.  23,  405-418  (1951). 

^Dunsiger,  A.  D.,  "High-Frequency  Acoustic  Echoes  Received  From  Simple 
Geometric  Shapes  with  Possible  Applications  to  Target  Recognition," 

J.  Sound  Vib.  13,  323-345  (1970). 
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target  echoes  were  a  sequence  of  pulses  whose  spacing  corresponded  to 
the  distance  between  edges  and  corners  of  the  target  as  viewed  from  the 
projector's  line  of  sight. 

Neubauer  and  Dragonette*  have  investigated  the  nature  of  the  wave- 
field  scattered  from  a  cylinder  at  broadside  incidence  by  means  of 
Schlieren  photography.  Their  work  has  contributed  significantly  to  our 
understanding  of  the  complex  interaction,  at  the  surface  of  the  target, 
between  compressional  waves  in  the  fluid  medium  and  compressional  and 
shear  waves  in  the  solid  target.  Flax®  has  shown  that,  because  of  these 
non-rigid  impedance  conditions,  the  responses  of  elastic  spheres  and 
cylinders  do  not  approach  a  limiting  value  at  high  frequencies  as  pre¬ 
dicted  by  a  physical  optics  approximation  for  the  scattering. 

Even  a  casual  familiarity  with  the  work  cited  above  suggests  that 
there  is  a  need  for  simultaneous  representation  of  impedance  effects  as 
well  as  aspect  and  frequency  dependence  of  a  target's  response.  Further 
more,  a  formalism  which  relates  boundary  impedance,  aspect  dependence 
and  frequency  dependence  with  experimental  geometry  (e.g.,  bistatic 
angle)  would  be  useful  in  interpreting  previous  and  future  measurements. 
In  the  author's  opinion,  this  formalism  need  not  be  mathematically  pre¬ 
cise  so  long  as  it  is  illuminating,  versatile  and  capable  of  numerical 
calculation.  It  is  in  this  vein  that  the  following  work  is  presented. 


^Neubauer,  W.  G.  and  Dragonette,  L.  R. ,  "Observation  of  Waves  Radiated 
from  Circular  Cylinders  Caused  by  an  Incident  Pulse,"  J.  Acoust.  Soc. 

Am.  48.  1135-1149  (1970). 

®Flax,  L.,  "High  ka  Scattering  of  Elastic  Cylinders  and  Spheres,"  J.  Acoust. 
Soc.  Am.  62,  1502-1503  (1977). 
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SECTION  II 

GENERAL  SCATTERING  THEORY 


Assume  that  a  target,  whose  surface,  A^^,  is  finite  and  closed,  is 
immersed  in  an  infinite,  homogeneous,  isotropic,  fluid  medium.  An 
acoustic  source,  P,  produces  a  compressional  wave  in  the  fluid,  and  a 
detector,  D,  senses  acoustic  waves  in  the  medium.  The  source-target- 
detector  geometry  is  shown  in  Figure  1.  An  arbitrary  origin  point,  0, 
has  been  chosen,  so  that  a  representative  point,  S,  on  the  target's 
surface  is  a  distance  r"  from  the  origin.  The  detector,  D,  is  located 
by  r  from  the  origin  point.  A  pseudo-origin,  T,  in  the  interior  of  the 
target  is  also  shown. 

All  transmitted  and  reflected  acoustic  waves  in  the  fluid  medium 
are  governed  by  the  wave  equation: 


v2H'(r,t)-l/c2  =  -  4np(r,t)  (1) 

where  p(r,t)  is  the  acoustic  source  distribution,  and  'f(r,t)  is  the 
acoustic  pressure  wave.  We  wish  to  know  the  form  of  ¥(r,t)  at  an  obser¬ 
vation  point,  r,  after  including  the  appropriate  scattering  effect  of 
the  target.  We  may  proceed  to  the  solution  by  the  method  of  Green's 
function.  Within  this  approach,  we  may  assume  an  impulse  source  and 
proceed  to  solve  for  the  time  dependence  of  H'Cr,!)  directly^®,  or  we  may 
assume  a  harmonic  time  dependence  and  solve  for  the  implicit  (and  equiva¬ 
lent)  frequency  dependence  via  the  Helmholtz  equation  (pp  803-811, 
Reference  10).  We  have  selected  the  latter  approach. 


^%orse,  P.  M.  and  Feshbach,  H.,  Methods  of  Theoretical  Physics,  McGraw 
Hill  Book  Company,  New  York,  pp.  834-848  (1953). 
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The  time  harmonic  assumption  is: 


'J'(r,t)  =  H'(r)e^'"*^  .  (2) 

Since  the  scattering  process  is  linear  for  infinitesimal  amplitude 
signals,  we  may  generate  the  time  series  for  any  transmitted  wave  shape 
by  Fourier  analysis  once  'l'(r)  has  been  determined. 


T(u))  'l'(r)e^‘"‘'du)  (3) 

where  T(u))  is  the  Fourier  transform  of  the  transmitted  signal.  With 
this  option  available,  we  now  return  to  the  wave  equation. 

Assume  a  harmonic  source  density  for  p(r,t)  similar  to  equation 
(2).  The  wave  equation  then  may  be  transformed  from  equation  (1)  to: 


V^4'(r)  +  k^'i'Cr)  =  -4n:p(r)  (4) 

where  the  time  dependence  has  been  eliminated.  Note  that  k  =  U)/c  where 
"c"  is  the  speed  of  the  acoustic  wave  in  the  fluid  medium.  Equation  (4) 
is  the  inhomogeneous  Helmholtz  equation.  By  definition,  a  Green's 
function  satisfies  the  Helmholtz  equation  for  a  unit  point  source  at  a 
location  r" .  Thus , 


V^G(r,r")  +  k^G(r,r")  =  -47t6(r-r")  (5) 

where  6(r-r")  is  the  Dirac  delta  function  and  G(r,r")  is  the  Green's 
function  relating  the  field  at  the  point  r  to  the  source  at  the  point 
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If  equation  (4)  is  multiplied  by  G(r,r")  and  equation  (5)  is  multi- 
plied  by  'i'Cr")  and  the  difference  taken,  the  result  is: 

G(r,r")V2'}'(r")  -  't'(>)V2G(r,>) 

=  4n[H'(>)6(r->)  -  G(r,>)p(>)].  (6) 

Note  that  the  arguments  in  equation  (4)  were  changed  from  r  to  r"  before 
multiplication.  Equation  (6)  is  qui.te  general  and  applies  for  any  points 
r  and  r"  in  a  common  volume.  We  now  enclose  the  volume  containing  all 
the  source  and  field  points  of  interest  by  two  surfaces  and  A2  as 
shown  in  Figure  1.  Surface  Aj^  is  the  target  surface  while  surface  A2  is 
at  a  large  but  arbitrary  distance  from  all  sources,  targets,  and  field 
points  of  interest.  Equation  (6)  may  then  be  integrated  over  the  volume 
"V"  bounded  by  Aj^  and  A2. 

^  S  (G(r,>)V2H'(^)-'l'(>)V2G(r,r")]dV(r") 

=  H'(r)  -  J  G(r,^)p(^)dV(^).  (7) 

V 

Note  that  the  integral  property  of  the  Dirac  delta  function  has  been 
used  to  produce  H'(r)  on  the  right  hand  side  (RHS)  of  equation  (7). 

The  volume  integral  on  the  left  hand  side  (LHS)  of  equation  (7)  can  be 
converted  to  a  surface  integral  via  Green's  Theorem.  Thus, 


J  [G(r,r")^(>)  -  ‘l'(r")^G(r,r")]-dl(> 
A  ■ 


+  5^  J  [G(r,>)^(>)  -  >Kr")^G(r,r")]-dl2(r" 


=  H'(r)  -  ;  G(r,>)p(?")dV(>) 
V 


(8) 
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where  di^j^(r")  and  are  shown  in  Figure  1.  The  contribution  of 

the  surface  integral  over  A2  can  be  made  vanishingly  small.  Various 
authors^^’^^  present  different  arguments  in  this  regard,  but  none  is 
mathematically  precise.  We  simply  shall  state  that  the  contribution 
from  surfaces  other  than  the  target  can' be  made  vanishingly  small  by 
careful  selection  of  the  experimental  conditions.  When  this  can  not  be 
arranged,  e.g.,  some  real  world  operational  situations,  the  result  is  a 
detectable  contribution  from  A2  known  as  simple  reverberation.  For  the 
purposes  of  this  description,  we  ignore  reverberation  so  that  equation  (8) 
becomes: 

^  S  [G(r,^)^(r■')  -  4'(^)^G(r,^)]•dX^(^) 


=  H'(r)  -  ;  G(r,>)p(>)dV(>)  .  (9) 

V 

The  volume  integral  on  the  RHS  of  equation  (9)  can  be  identified  by 
assuming  the  target  is  absent.  Hence,  the  integral  over  A^  is  zero  and: 

‘H(r)  =  f  G(t,P')p(r")dV(r")  .  (10) 

V 


The  equation  simply  states  that  the  wave  detected  at  r  originates  solely 

from  the  acoustic  sources  at  r".  Hence,  the  integral  on  the  RHS  of 

equation  (10)  is  the  incident  wave.  Equation  (9)  may  be  rewritten  in 

terms  of  the  incident  wave,  ib-  (r),  and  the  scattered  wave,  ib  (r)  as: 

’  ^inc^ 

'}'(r)  =  H'.  (r)  +  ¥  (r)  (11) 

'  me'  s'  ' 

where 

(r)  =  ;  G(r,r")p(>)dV(^)  (12) 

me  y 


^^Marion,  J.  B. ,  Classical  Electromagnetic  Radiation,  Academic  Press, 
New  York,  1965. 

^^Jackson,  J.  D.,  Classical  Electrodynamics,  John  Wiley  &  Sons,  Inc., 
New  York,  1962. 
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'-4i 


and 

't's(r)=^J  [G(r,>)^(>)-4'(^)^G(r,^)]-d^^(>)  .  (13) 

In  order  to  agree  with  the  normal  convention  for  outward  pointing  sur¬ 
face  normals,  the  direction  of  d^j^(r")  will  be  reversed  (accompanied  by 
the  appropriate  sign  change)  so  that  it  points  outward  from  the  target. 

For  purposes  of  simplification,  the  subscripts  will  be  dropped  since  it 
will  be  understood  henceforth  that  we  are  describing  the  scattered 
wavefield,  4*g(i^)>  arising  from  the  interaction  of  the  incident  wave  with 
the  target  surface  only.  Hence, 


/  [H'(^)^G(r,r")  -  G(r,r'')^(^)]'d^(r")  .  (14) 

A 

Equation  (14)  is  the  Kirchhoff  diffraction  equation.  It  is  quite 
general  and  should  not  be  confused  with  the  "Kirchhoff  approximation" 
which  is  concerned  with  an  approximate  solution  to  equation  (14)  involving 
a  number  of  assumptions. 
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SECTION  III 


EMPIRICAL  SCATTERING  MODEL 


If  the  value  of  the  wavefield  and  its  gradient  were  known  at  each 
point  of  the  target's  surface,  then  equation  (14)  could  be  used  with  a 
spherical  wave  Green's  function  to  solve  for  the  scattered  wavefield  at 
any  desired  point  r.  Since  we  do  not  know  the  value  of  the  wavefield, 
let  us  assume  that  it  can  be  represented  as  a  linear  combination  of 
component  waves;  e.g.,  the  incident  wave,  diffracted  waves,  and  multiple 
scattered  waves. 


'J'Cr")  =  I 

j=l  J  J 


where  the  fljCr")  are  complex  distribution  functions  specifying  the 
amplitude  and  phase  of  each  of  the  "j"  cotoponents  at  each  point  on  the 


target's  surface. 


W(r")  =  I  [a.(r")^.(r")  + 'l'.(r")^a.(>)]  . 

j— J  J  J  J 

Substituting  equations  (15)  and  (16)  into  equation  (14)  yields: 

H'(r)  =  Z  4  [aj(>)H'j(?")^G(r,>)  -  aj(?")G(r,r")^j(^) 
J  ^ 


G ( r ,  r"  )¥ Cr "  )^a .  (r" )  ]  •  d^ ( r " ) } 
J 
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Note,  however,  that  each  C(j(r")  is  a  surface  function,  varying  only 
parallel  to  the  surface  at  any  point.  Hence,  ^0(.(r")  is  parallel  to  the 


surface  and  perpendicular  everywhere  to 


dUP).  ^ 


Thus , 


(x") ‘d^Cr")  =  0  for  all  "j"  . 
Therefore,  we  rewrite  equation  (17)  as: 

->  » 

H'(r)  =  I  Hf-Cr) 

j=l  J 


where 


J  J  J 


r,r'')-G(r,r")^j  (r'')]  ‘d^Cr")  . 


(18) 


(19) 


(20) 


The  task  which  remains  is  to  postulate  the  nature  and  origin  of  the 

component  waves.,  H'.(r"),  and  the  distribution  functions,  C(.(r"),  through 
•3  3 

physical  reasoning  and/or  experimental  measurements.  As  each  component 

is  added  to  the  summation,  its  contribution  can  be  compared  to  experi¬ 
mental  measurements  to  determine  its  significance. 

Before  proceeding  further,  we  wish  to  add  an  additional  simplifica¬ 
tion.  We  will  assume  that  the  target  is  in  the  farfield  of  both  the 
source,  P,  and  the  detector  at  D  in  Figure  1. 

Now, 


....  -iki?-^. 

G(r,r")  = 


-r" 


r-r 


(21) 


SD’ 


Using  the  farfield  approximation  and  noting  that  rg^  =  (see 


Figure  1) ,  then 


G(t,E")  ^'^^*G(e')  =  i- — - )  e 

’^TD 


(22) 
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where  k2  is  a  wave  vector  parallel  to  as  shown  in  Figure  1.  Care 
now  must  be  exercised  with  respect  to  the  gradient  operators  in  equation 

(20) .  The  gradient,  ^G(r,r"),  is  with  respect  to  r"  while  VG(r')  is 

with  respect  to  r' .  Noting  the  reversal  of  signs  of  r"  and  r*  in  equations 

(21)  and  (22),  we  see  that: 


We  also  will  assume 


-ik|r"-rj  I  -iklr.gl 


i-> 


|r  ,..g 

and  use  the  farfield  approximation 


-il-r-T 

H'.tr")^-^-  4'|(r')  =  - ~)  e  ^ 

J  J  r.^ 


(23) 


(24) 


(25) 


where  is  a  wave  vector  from  the  point  "j"  to  the  target  origin  point 
"T"  parallel  to- the  vector  r.™.  In  Figure.!,  this  is  illustrated  for 
the  incident  wave  from  the  point  "P"  where  kj  =  and  r^j  =  rp^. 

Noting  the  similarity  in  sign  for  r"  in  equation  (24)  and  r'  in 
equation  (25),  we  see  that 


^.(r")  =  ^.(r')  .  (26) 

J  <J 

Now  from  equations  (25)  and  (22) , 

^.(^')  =  -il  'F.(r')  (27) 

J  J  J 

and 

^G(r')  =  i^2G(r’)  •  (28) 
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Substituting  equations  (22),  (25),  (27),  and  (28)  into  equation  (20) 
yields; 


y.(?) 


-ik(r^jj+r.^) 


i(^2"^j)*r' 


]  •  dX(?'). 


(29) 


Let  us  define  $2j  ~  ^2~^j  =  ndA(r')  where  n  is  an  outward 

pointing  unit  normal. 

Then, 


4'j(r)  = 


■ik(r,jj,jj+rj.j,) 


4n 


-)1  ;  art') 

A  ■’  l''2jJ 


)  e 


i^  -r' 

mP)  . 

(30) 


The  factor  in  brackets  contains  the  experimental  information  regarding 
the  distances  from  the  incident  wave  source  "j"  to  the  target  origin  and 
from  the  target  origin  to  the  point  of  observation  "D".  The  integral, 
on  the  other  hand,  appears  characteristic  of  the  target  and  its  relative 
orientation  with  respect  to  source  and  detector  directions. 


Before  examining  the  integral  more  closely,  the  farfield  assumption 
merits  some  discussion.  For  sources  external  to  the  target,  the 
farfield  assumption  is  not  very  restrictive.  For  waves  which  may  arise 
from  multiple  reflections  in  or  on  the  target,  this  assumption  appears 
invalid  in  principle.  However,  because  we  are  at  liberty  to  adjust  the 
aj(r')  distribution  functions  in  amplitude  and  phase,  we  have  sacrificed 
no  generality  in  principle.  We  have  made  no  real  assumption  and  hence 
produced  no  apparent  simplification.  The  simplification  is  in  form  only 
as  manifested  by  the  form  of  the  integral  in  equation  (30) . 


Let  us  define 
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FA(ty)  = 


dA(r' ) 


(31) 


where  structure  factor  of  target  "A".  Note  that  the 

structure  factor  and  the  factor  inside  the  braces  of  the  integral  are  a 
Fourier  transform  pair.  The  factor  in  braces  may  be  called  the  "j"th 
component  of  the  target  characteristic  function.  Because  of  the  Fourier 
relationship,  we  may  construct  a  wave  space  for  which  is  conju¬ 

gate  to  the  real  space  which  is  the  domain  of  the  characteristic  function. 
This  concept  originated  over  50  years  ago  in  x-ray  diffraction^^’ and 
has  been  discussed  more  recently  by  Lewis^®  in  connection  with  physical 
optics  inverse  diffraction.  The  Ewald  construction,  which  relates  the 
incident  wave  vector  kj  =  kj^  and  the  scattered  wave  vector  k2  to  the 
argument  ^2]^  of  ^A^^21^’  discussed  briefly^®  but  will  be 

included  for  completeness  in  the  description  of  wave  space  in  the  next 
section. 


^^Ewald,  P.  P.,  "Zur  Theorie  der  Interferenzen  der  Rontgenstrahlen  in 
Kristallen,"  Physik.  Z.  14,  465  (1913). 

^^Pepinsky,  R.  and  Vand,  V.,  "Crystallography  and  X-Ray  Diffraction," 
Handbook  of  Physics,  E.  V.  Condon  and  H.  Odishaw  (eds.),  McGraw-Hill 
Book  Company,  Inc.,  New  York,  pp.  8-1  to  8-23  (1958). 

^®Lewis,  R.  M.,  "Physical  Optics  Inverse  Diffraction,"  IEEE  Transactions 
on  Antennas  and  Propagation  17,  308-314  (1969). 

^^Nelander,  J.  C.,  "Modeling  Acoustic  Scattering  from  Elastic  Cylinders," 
NCSC  Technical  Note  477,  March  1977. 


14 


15 


NCSC  TM-275-79 


and 


+  kyY  +  k^2  .  (37) 

The  basis  vectors  in  either  wave  or  real  space  need  not  be  orthononnal^^ 
and  usually  are  not  when  symnietry  conditions  are  more  important  than  the 
simplification  produced  by  orthonormality.  For  our  purposes,  we  shall 
assume  orthonormality. 

The  conjugate  relationship  between  a  real  space  characteristic 
function  and  a  wave  space  structure  factor  is  illustrated  in  Fi.iure  2. 

The  Fourier  relationship  between  the  two,  equation  (32),  ensures  that 
complete  knowledge  of  one  is  equivalent  to  complete  knowledge  of  the 
other.  Equation  (30)  indicates  that  we  can  gain  knowledge  about  the 
structure  factor  (and,  hence,  the  target  characteristic  function)  by 


Characteristic  Function  ^  Structure  Factor  F^  (k) 

Ikl 


FIGURE  2.  REAL  SPACE  -  WAVE  SPACE  RELATIONSHIP 


^^ibid 
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measuring  the  acoustic  wavefield  scattered  by  a  target.  Several  impor¬ 
tant  points  warrant  further  discussion.  For  example:  (1)  how  is  the 
characteristic  function  related  to  such  familiar  target  parameters  as 
shape,  size,  and  orientation?  (2)  what  parts  of  the  structure  factor 
can  be  measured  in  a  given  source-target-detector  geometry  utilizing  a 
specific  transmitted  waveform?  (3)  is  the  structure  factor,  and  hence 
target  characteristic  function,  unique;  i.e.,  independent  of  source- 
target-detector  geometry?  (4)  is  there  a  minimal  set  of  measurementr- 
which  can  be  made  which  will  identify  a  given  class  of  targets  with  a 
high  probability  of  success?  The  answer  to  the  fourth  question  can  be 
attempted  only  after  completing  a  modeling  and  measurement  effort. 
Information  which  will  be  useful  in  answering  the  other  questions  is 
provided  here. 

EWALD  CONSTRUCTION 

The  Ewald  construction^^  can  be  used  to  relate  experimental  param- 
eters  to  the  range  of  F(k)  values  which  are  observable  by  that  experiment. 
Assume  that  the  incident  plane  wave  is  in  the  direction  In  wave 
space,  a  sphere  of  radius  |  =  2n/X  can  be  drawn  such  that  its  diameter, 
along  uhe  direction,  passes  through  the  origin  of  wave  space  coordinates 
so  that  the  sphere  is  tangent  to  the  origin.  The  scattered  wave  vector, 

-J* 

is  drawn  from  the  center  of  this  sphere  in  the  scattered  wave  direction. 
The  resulting  constructon  is  shown  in  Figure  3(A).  Note  that  the  vector 
fc  =  k2-kj^  has  its  terminus  on  the  sphere  regardless  of  the  direction  of 
1^2*  Hence,  for  a  given  target  orientation  with  respect  to  only 
those  values  of  F(^)  lying  on  the  surface  of  the  sphere  are  observable 
regardless  of  the  position  of  the  detector.  If  a  transmitted  signal  has 
a  finite  bandwidth  (rather  than  monochromatic),  then  two  spheres  with 
radii  1^,  |(max)  and  11^, | (min)  may  be  drawn  as  shown  in  Figure  3(B)  .  The 
high  frequency  limit  of  the  incident  wave,  |k^|(max),  defines  the  larger 
sphere  while  the  low  frequency  limit,  |kj^|(min),  specifies  the  smaller 


^^ibid 
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(A)  Single  Frequency  (B)  Broadband 

FIGURE  3.  EWALD  CONSTRUCTION 


sphere.  Observable  values  of  F(k)  are  contained  in  the  volume  inside 
the  larger  sphere  but  outside  the  smaller  sphere.  The  line  of  F(^) 
values  which  can  be  observed  with  a  broadband  incident  wave  and  a  back- 
scattering  configuration,  i.e.,  i.2  ~  is  shown  in  Figure  3(B).  With 
the  help  of  Figure  3(B)  and  a  little  imagination,  one  can  see  that  the 
same  F(1c)  value  can  be  observed  by  more  than  one  experimental  configura¬ 
tion.  But  these  values  will  occur  for  different  frequencies.  It  should 
be  emphasized  that  these  would  be  the  same  F(^)  values,  not  just  equal 
F(1c)  values.  This  occurs  because  of  the  multitude  of  ways  that  the  same 
vector  k  can  be  constructed  from  different  kj^  and  k2  vectors.  Note  that 

^1  ’  determines  the  frequency  at  which 
the  F(Ic)  value  will  be  observed. 
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Although  the  above  relationships  are  literally  true  for  diffraction 
where  equation  (32)  represents  the  total  structure  factor,  the  relation¬ 
ships  for  acoustic  scattering,  where  multiple-scattered  waves  may  be  of 
considerable  significance,  are  bound  to  be  more  complex.  Suffice  it  to 
say  that  caution  should  be  exercised  when  comparing  backscattering 
measurements  to  bistatic  measurements  for  the  same  target.  Some  polar 
plots,  such  as  those  by  Barnard  and  McKinney^,  which  display  the  scatter¬ 
ing  cross  section,  |F(^)| where  the  projector  and  target  are 

stationary  while  the  detector  is  moved  in  a  circle  around  the  target, 
are  variable  bistatic  angle  measurements.  The  F(k)  values  so  measured 
would  fall  on  a  circle  on  the  surface  of  a  single  sphere  such  as  that 
shown  in  Figure  3A.  If  the  same  target  were  examined  by  backscattering 
measurements,  using  the  same  frequency,  the  series  of  F(k)  values  so 
measured  would  trace  out  a  circle  of  radius  2|kj^|whose  center  is  the 
origin  of  the  wave  space  coordinate  system  rather  than  the  center  of  the 
Ewald  sphere  in  Figure  3(A).  Whether  or  not  there  exists  an  operator 
such  that: 

%  ^^^^^bistatic^  "  ^backscattering 


for  acoustic  scattering  is  not  known.  It  is  the  author's  opinion  that  a 
thorough  understanding  < 
of  acoustic  scattering. 


thorough  understanding  of  is  equivalent  to  a  thorough  understanding 


INCLINATION  FACTOR 


One  final  point  is  worth  amplification.  The  factor. 


•  n 


in  equation  (32)  is  simply  a  variant  of  the  Stokes  ir  'lination  factor  of 
physical  optics.  Consider  the  reflection  of  an  incident  wave  k^^  from  a 
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planar  surface,  whose  unit  normal  is  h,  into  the  specular  direction  k 


^1  = 


and  =  6^,  it  is  easy  to  see 


as  shown  in  Figure  4.  Since 
that  k  =  k  “k,  is  parallel  to  h  and  hence  the  inclination  factor  as 

S  X 

maximum;  i.e., 


t 


1 


(3S) 


For  directions  other  than  specular,  i.e.,  k,  k  ,  the  inclination 

X  S 

factor  specifies  the  appropriate  weight  for  the  reflectivity  of  that 
point  into  the  direction  specified  by  inclination  factor  is 

non-negative  and  is  minimum  when  k^  is  parallel  to  the  surface  and  ^2  is 
parallel  or  anti-parallel  to  When  applying  the  inclination  factor, 
one  need  not  be  concerned  with  deviations  from  the  specular  direction 
but  only  with  the  relationship  between  ^  and  h. 


(Text  Continued  on  Page  22) 
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SECTION  V 

RIGID  CYLINDER  SCATTERING  -  FIRST  LEVEL  APPROXIMATION 

GENERAL  APPROACH 

In  order  to  illustrate  the  utility  of  the  empirical  model 
derived  in  Secton  III,  let  us  consider  scattering  from  a  rigid,  finite 
cylinder  at  an  arbitrary  orientation  or  aspect  angle.  Objects  with 
edges,  such  as  the  finite  cylinder,  pose  particularly  difficult  problems 
in  scattering  theory  because:  (1)  the  edges  produce  multiple  scattering 
effects,  and  (2)  the  edges  are  the  loci  of  discontinuities  in  geometry 
and  hence  make  integration  of  wavefields  across  these  areas  difficult. 

A  simple  scattering  model  by  Freedman^’  cites  these  discontinuities 
as  the  principle  source  of  scattering.  Freedman's  image  point  model  of 
acoustic  scattering  is  a  scalar  model  in  which  the  target's  surface  is 
viewed  from  a  common  source-detector  position  (backs cattering) .  Dunsiger^ 
has  noted  that  the  integral  of  the  image  point  model  is  difficult  to 
evaluate  because  of  discontinuities  in  the  integrand.  These  discontinu¬ 
ities  are  the  image  points.  Although  the  structure  factor  for  our 
empirical  model  given  in  equation  (31)  is  similar  to  Freedman's  integral, a 
judicious  choice  of  integration  procedure  can  avoid  mathematical  difficulties 
Note  also  that  equation  (31)  is  a  vector  rather  than  scalar  representation. 
The  integration  now  will  be  illustrated  for  the  rigid,  right  circular 
cylinder  shown  in  Figure  5. 

The  radius  of  the  cylinder  is  "a,"  and  its  length  is  "L."  The 
cylinder  axis  is  aligned  parallel  to  the  real  space  z-axis.  We  divide 
the  surface  of  the  cylinder  into  two  parts:  Sj^,  the  end  disc,  and  S2, 


^^Freedman,  A.,  "A  Mechanism  of  Acoustic  Echo  Formation,"  Acustica  12, 
10-21  (1962). 

(Text  Continued  on  Page  24) 
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Y 


FIGURE  5.  GEOMETRY  FOR  SCATTERING  FROM  A  RIGID 
RIGHT  CIRCULAR  CYLINDER 
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the  curved  surface.  We  will  consider  only  the  first  level  approxi¬ 
mation;  i.e.,  only  the  incident  wave  from  the  source  is  considered  and 
all  body-diffracted  and  multiple-scattered  waves  are  ignored. 

The  structure  factor  is: 


FCi!)  =  X  «(?)  fel  dA(?) 

Ui 

The  distribution  function  is: 

a(r)  =  I  0[-k, -n  ]  6[k.,‘n  ] 

'■  1  m-*  ■■  2 

m=l 

where  0(x)  =  1,  x>0;  0(x)  =  0  x<0 


(AO) 


(Al) 

(A2) 


and  "ra"  denotes  the  various  analytic  surfaces  of  the  target.  Hence, 
equation  (Al)  states  that  only  those  surfaces  observable  from  both  the 
source  and  detector  are  of  significance.  We  may  confine  our  investiga¬ 
tion  to  the  x-z  quadrant  of  Figure  5,  because  of  symmetry.  Further,  we 
require  fc2  to  be  coplanar  with  k^^  and  the  cylinder  axis.  This  imposes 
some  restriction  on  the  generality  of  our  solution.  Equation  (Al)  is 
nonzero  for  and  one  half  the  curved  side  surface  (S2)  of  the  cylinder. 
We  may  write: 


F(]|)  =  F„  (S)  +  F.  (^) 
^1  ^2 


(A3) 


where  F„  and  F^  are  the  structure  factors  of  the  end  disc  and 

hemi-cylinder  respectively.  Each  structure  factor  may  be  evaluated 
separately  and  summed  to  produce  the  total  F(5). 


CYLINDER  END  DISC  STRUCTURE  FACTOR 


For  the  end  disc: 
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Sj  IJI  1 


Now  from  Figure  6,  we  see  that: 

r  =  2  2  +  X 

and,  dA(r)  =  2ydx  =  2(a^-x^)^  dx  . 
Thus, 


.  +a  j  ,  ^  i(k  x+k  z) 

Fo  tk)  =  /  ;  6(2-^)  2(a2-x2)^  (k^)  e  ^ 

h  z  -a  2  1^1 


dxdz 


(43) 


(44) 

(45) 


(46) 


Note,  from  Figure  5,  ^‘n  =  ^‘z  =  k^  for  the  surface  Sj^. 

Thus  equation  (46)  can  be  integrated  over  dz  and  simplified  to: 


F(^) 


+a 

J 


-a 


L  ik  X 

(a2-x2)^  e  dx. 


The  integral  in  equation  (47)  can  be  simplified  as  shown  below. 


+a  1  ik  X  a  , 

J  (a2-x2)'2e  ^  dx  =  2  X  (a2-x2)'2cos(k^x)dx. 


(47) 


(48) 


Let,  t=-;  dt=— ;  x=0  t=0;  x=a  ^  t=l. 

d  d 

Thus  the  integral  becomes: 

1  , 

2a2  X  (l“t2)\os[(k  a)t]dt. 

0  * 

The  integral  (49)  is  a  form  of  the  first  order  Bessel  function.^® 
Thus, 


(49) 


^’Abramowitz,  M.  and  Stegun,  I.  A.,  Handbook  of  Mathematical  Functions, 
Dover  Publication,  Inc.,  New  York,  1965. 
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«  ik  I,/2  k  JiCk  a) 

F  (k)  =  (2na=‘e  )  f-Vl-l 

1  {k|  X 

where  Jj^(k^a)  is  a  first-order  Bessel  function, 

-  (■>-)  J„  ITT^) 


and  "F"  is  the  gaouna  function. 


CYLINDER  SIDE  SURFACE  STRUCTURE  FACTOR 

In  order  to  evaluate  the  structure  factor  for  the  hemi-cylinder, 
consider  the  diagram  in  Figure  7.  Note  that  r  =  p+z  and  p  =  x+y.  Now 
the  unit  normal  to  the  surface  is: 


A  ->  ->  •> 

n  =  fi  = 
a  a 


Thus , 


^  k  X 

k»n  _  k*x  +  k*y  _  x 
k|  |k|  |k| 


(51) 


since  k*y  =  0  because  k  lies  in  the  x-z  plane  (see  above).  The  dif- 


or 


area  element  in  Figure  7  can  be  written  as: 

dA  =  ad(|idz 

adxdz 

(for  a>x>0) 

(52a) 

(a^-x^)^ 

_  1  adxdz 

(for  0<x<a)  . 

(52b) 

(a2-x2)^ 

The  two  forms  of  equation  (52)  are  needed  for  the  positive  and  negative 
values  of  "y."  Thus,  we  may  write: 
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1  ,  ik  at 

=  1  +  ik  a  J  (1-t^)^  e  ^  dt. 
^  0 


In  order  to  identify  the  integral  in  equation  (58),  we  write  it  in  the 
following  form: 


f  (1-t^)^  e  ^  dt  =  /  (1-t^)^  cos((k  a)t]dt 


+  i  /  (1-t^)^  sin  [(k  a)t]dt  . 

0  ^ 


The  first  integral  on  the  RHS  of  equation  (59)  is  a  form  of  the  first- 
order  Bessel  function^®  (see  also  equation  (49))  while  the  second  inte¬ 
gral  on  the  RHS  of  equation  (59)  is  a  form  of  the  first-order  Struve 
function^®,  Hj^(z),  where: 


n  /l,,^2n 


H  (z)  =  (^Z)^  Z  (~1)  (^Z) 

H^tzj  1  r(n+3/2)  r(n+ 


where  F  is  the  gamma  function. 

Equation  (58)  now  can  be  put  in  the  form: 


=  [l-^^(k^a)]  +  i[|j3^(k^a)] 


Using  equations  (56)  and  (61)  in  equation  (55)  we  find: 

2  |kl  (_z_) 


^®ibid 
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COMBINED  STRUCTURE  FACTOR 

la  order  to  combine  F<;  (k)  with  F„  (k) ,  we  rewrite  the  exponential 
factor  in  equation  (50) , 


e  =  cos  (-|-)  +  i  sin  (-|-)  .  (63) 

Equation  (50)  then  becomes: 

Fg  it)  =  27Ca^(4f)[^^-][cos  i^)  +  i  sin  i^)]  (64) 

Combining  equations  (62)  and  (64)  yields: 

F(^)  =  F.  it)  +  F.  it) 


2na2k,  k,L  2naLk^  k,L  „ 

(—•cos  (^)Ji(k,a)  .  ^  sin  (-) [j-HiCk/) ] } 


2na^k  k  L  27taLk  k  L 

+  i{^;^ - -  sin  (-f-)Ji(k  a)  +  — - -  sin  (-|-)  J-,(k  a)}  . 

Ikjk^a  2  ^  *  ll^lk^L  2  1  X  (gg) 


Thus,  the  total  structure  factor,  to  this  first  level  approximation,  for 
the  rigid,  right  circular  cylinder  is: 


27ta  z 

F(k)  -  -r:  iiT  cos(-^ 


V„2 


jfj  fk"  ir  <— 


k  k  k  L 

+  i  sin(-|-)  Jj^(k^a)} 

X  z 


where  J,  (k  a)  is  first-order  Bessel  function  and  H-,  (k  a)  is  a  first- 

JL  X  X  X 

order  Struve  function  defined  in  series  form  by  equations  (50a)  and 
(60),  respectively. 
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SECTION  VI 

RIGID  CYLINDER  SCATTERING  CROSS  SECTION 


The  structure  factor,  equation  (66),  is  valid  for  both  bistatic 
scattering  and  backscattering,  to  the  level  of  approximation  indicated 
by  equations  (41)  and  (42),  as  long  as  kj^,  and  the  cylinder  axis  are 
coplanar.  The  nature  of  the  structure  factor  can  be  examined  by  plotting 
the  dimensionless  scattering  cross-section, 

0(D  =  [C-F(0]2  (67) 

where  -  |^|a.  The  dimensionless  parameter  ^  is  proportional  to  the 
ratio  of  the  nominal  target  dimension  (in  this  case,  the  cylinder  radius 
"a")  to  the  wavelength  of  the  acoustic  waves  used  to  provide  the  scattering 
information.  Hence,  a  particular  value  of  0(0  is  representative  of  a 
large  target  at  low  frequency  or  a  small  target  at  high  frequency  when 
similar  scaling  factors  are  applied  to  the  target  dimension  and  wavelength. 

A  three-dimensional  plot  in  wave  space  of  the  scattering  cross-section 
for  a  rigid  cylinder  whose  radius  to  length  ratio  is  6  is  shown  in 
Figure  8.  The  ordinate  ij  Figure  8,  Q(0 ,  is  measured  in  decibels 
relative  to  the  maximum  value  in  the  plot.  The  axis  corresponds  to 
the  end-on  aspect  of  the  cylinder  while  the  axis  is  the  broadside 
response.  Only  one  quadrant  of  the  function  fl(0  is  shown  since  the 
others  are  similar  by  reason  of  symmetry.  A  shaded  contour  plot  of  Q(0 
showing  all  four  quadrants  is  shown  in  Figure  9.  Generally  speaking, 
the  darker  the  region  in  Figure  9,  the  larger  the  f2(t)  value.  However, 
caution  should  be  exercised  when  interpreting  Figure  9  since  it  is  a 
black  and  white  reproduction  of  a  color-coded  display. 

(Text  Continued  on  Page  35) 
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FIGURE  9.  SHADED  CONTOUR  PLOT  SHOWING  THE  FOUR  QUADRANTS  OF 
THE  RIGID  CYLINDER  SCATTERING  CROSS-SECTION 
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When  interpreting  either  Figure  8  or  Figure  9,  one  should  imagine 
that  the  target  is  located  at  the  origin  of  coordinates  of  each  plot 
with  the  cylinder  axis  oriented  along  the  axis.  The  response  which 
would  be  observed  in  a  given  backscattering  situation  can  be  found  by 
drawing  a  radial  line  from  the  origin  of  the  plot  in  the  direction 
corresponding  to  the  target  viewing  or  aspect  angle.  The  appropriate 
response  is  located  on  this  line  at  a  distance  from  the  origin  cor¬ 
responding  to  An  times  the  ratio  of  the  radius  "a"  to  the  wavelength 
"X."  Thus,  the  responses  for  higher  frequencies  are  located  at  larger 
radial  distances  from  the  origin.  The  maximum  values  of  ^  =  47t(a/X) 
shown  in  Figures  8  and  9  is  25. 

A  few  comments  on  the  interpretation  of  Figure  8  are  worthwhile. 
First,  the  scattering  cross-section  increases  monotonically  with  frequency 
for  constant  values  of  the  structure  factor,  F(t).  In  a  linear  plot, 
this  increase  would  be  quadratic.  It  arises  because  of  the  factor  in 

•  This  increase  can  be  seen  along  the  t  axis  of  Figure  8  where 
F(0  is  constant.  The  factor  is  a  part  of  the  definition  of  scattering 
cross-section  since  0(0  conventionally  is  defined  in  terms  of  energy 
scattered  into  a  given  solid  angle.  In  order  to  compensate  for  the 
spherical  spreading  of  the  energy  of  a  scattered  wavefield  (remember, 
since  'l'(r)  is  proportional  to  F(k),  or  in  terms  of  our  dimensionless 
parameter  F(0>  the  energy  density  F^(0  must  be  multiplied  by  a  factor 
proportional  to  the  area  of  a  sphere  in  wavespace,  i.e.,  the  radius 
squared,  0)* 

Second,  the  ripple  contours  whose  variation  is  parallel  to  the 
axis  (i.e.,  those  peaks  and  troughs  which  are  parallel  to  arise  pre¬ 
dominately  from  the  end  disc  of  the  cylinder.  The  spacing  of  the  peaks 
and  troughs  in  the  direction  is  an  inverse  function  of  the  diameter 
of  the  end  disc.  '  The  larger  the  diameter,  the  smaller  the  spacing. 

This  can  be  compared  to  the  more  rapid  fluctuations  or  smaller  spacing 
between  the  peaks  and  troughs  along  the  ^  direction  which  arise  from 
the  length  of  the  cylinder. 
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Third,  the  scattering  cross-section  and  its  slope  are  zero  at  the 
origin.  This  is  more  easily  seen  in  Figure  10  which  corresponds  to 
Figure  8  from  a  different  viewing  angle.  Figure  10a  is  a  logarithmic 
plot  (in  dB)  while  Figure  10b  is  a  linear  plot  of  the  region  neat  the 
origin  which  shows  the  zero  value  of  the  scattering  cross-section  and 
its  slope.  The  region  nearest  the  origin  in  wave  space  is  known  as  the 
Rayleigh  region.  The  scattering  cross-section  in  this  region  increases 
approximately  as  the  fourth  power  of  the  argument  The  rate  of  increase 
is  related  to  the  volume  of  the  target.  In  our  empirical  model,  it 
represents  the  "first  moment"  of  the  inclination  factor.  Consider  the 
general  expression  for  the  structure  factor,  equation  (40),  for  small 
values  of  the  argument  k. 


F(^->0)  =  Ja(r) 


fe]  dA(^) 

|k| 


Thus 


(68) 


(69) 


For  0((r)  =  1,  the  first  integral  in  equation  (69)  is  a  constant  equal  to 
the  average  value  of  the  inclination  factor.  The  second  integral  is 
proportional  to  k  ,  and,  because  of  the  r  factor  in  k*r,  represents 
what  might  be  called  the  first  moment  of  the  inclination  factor.  There¬ 
fore,  equation  (69)  can  be  written  as: 


F(^^O)  =  Cq  +  C^|]^ 


(70) 
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FIGURE  10b.  LINEAR  PLOT  SHOWING  THE  RAYLEIGH  REGION  OF 
THE  RIGID  CYLINDER  SCATTERING  CROSS-SECTION 
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where  Cq  is  the  zeroth  moment  or  average  value  of  the  inclination  factor 
and  is  the  first  moment.  It  can  be  argued  that  Cj^»Cq  since  the 
moment  arms  are  roughly  equal  to  the  average  target-center  to  target- 
surface  distance  (compare,  for  example,  the  area  of  a  spherical  shell  to 
the  first  moment  of  the  area  of  a  spherical  shell).  Thus, 

[F(]|->0)]2  =  eg  +  2CqCJ^|  +  Cl\t\^  (71) 

or  in  terms  of  the  dimensionless  parameter 
[F(^0)]2  =  c2  +  2CQC^t  +  . 

Therefore, 

fi(^O)  =  . 

From  equation  (73)  we  can  see  chat  the  amplitude  and  slope  of  fi(0  at 
the  origin  are  zero.  For  small  values  of  the  argument  the  behavior 
of  fi(0  goes  approximately  as  the  fourth  power  of  t  because  of  the 
dominance  of  over  Cq. 

For  cases  where  a(r)  f  1,  our  analysis  obviously  must  be  modified. 

For  example,  if  a(r)  can  be  approximated  by  a  power  series  in  "r,"  then 
the  constants  in  equation  (70)  are  the  zeroth  and  first  moments  of  the 
factor  composed  of  the  product  a(r)  and  the  inclination  factor. 

Returning  to  Figure  8,  we  note  that  except  for  the  periodicity  of 
the  crests  and  troughs  in  the  ^  and  ^  directions  and  the  behavior  of 
fi(P  near  the  origin,  there  are  few,  if  any,  features  of  interest  in  the 
scattering  cross  section  of  the  rigid,  finite  cylinder.  This  raises  a 
subtle  but  important  point.  The  plots.  Figures  8  and  9,  of  the  scat¬ 
tering  cross  section  permit  one,  at  a  single  glance,  to  note  any  signifi¬ 
cant  features  which  might  be  of  use  in  target  classification.  If  features 
are  present,  their  variation  as  a  function  of  frequency  and  target 


(72) 

(73) 


39 


NCSC  TM-275-75 


aspect  can  be  noted  easily.  The  probability  that  significant  features 
are  present  for  real  targets  appears  high  based  upon  results  published 
in  the  literature.®’®  As  we  will  see  in  Section  VIII,  incorporating 
impedance  boundary  conditions  into  our  model  is  sufficient  to  produce 
interesting  and  significant  features  in  the  scattering  cross-section  of 
the  finite  cylinder.  A  collection  of  plots  similar  to  Figure  8,  but  for 
a  class  of  targets  of  interest,  could  be  quite  useful  in  the  design  of  a 
sonar  whose  parameters  (such  as  bandwidth  and  center  frequency)  are 
selected  in  order  to  optimize  target  classification  or  detection. 


®ibid 

®ibid 
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SECTION  VII 

IMPEDANCE  BOUNDARY  CONDITION 


GENERAL  CONSIDERATIONS 

In  order  to  develop  a  second  level  approximation  to  the  model, 
one  which  incorporates  the  acoustic  impedance  properties  of  the  target 
surface,  let  us  consider  the  boundary  conditions.  Brekhovskikh^®  has 
derived  reflection  and  refraction  coefficients  for  the  case  of  a  compres- 
sional  plane  wave  in  a  fluid  medium  incident  upon  an  infinite,  plane, 
solid  surface  by  solving  the  boundary  value  problem.  The  geometry  and 
parameters  are  shown  in  Figure  11.  The  densities  and  wave  speeds  in  the 
two  media  are  indicated  in  the  diagram.  As  Brekhovskikh  shows,  the 
incident,  reflected,  and  refracted  wave  vectors  all  lie  in  a  common 
plane  (the  X-Z  plane  in  Figure  11)  as  a  consequence  of  the  boundary 
conditions.  For  an  incident  wave  of  unit  amplitude,  the  complex  ampli¬ 
tude  of  the  wave  reflected  back  into  the  fluid  medium  is : 


(74) 


(75) 

where  the  parameters  are  specified  in  Figure  11. 


where 


Zjj  cos2(20g)  +  Zj,  sin2(28g)  -  Zg 


Zg  cos^(20j,)  +  Zg  sin^(29g)  +  Zg 


•  Z  = 

’  ^D  cos  0 


_  .  z  - 

D  ’  E  cos  0g 


^^Brekhovskikh,  L.  M., 
1960. 


Waves  in  layered  Media,  Academic  Press,  New  York, 


(Text  Continued  on  Page  43) 
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FLUID 

Pp  -  Density 

Cp  -  Compressional 
Wave  Speed 


\  / 


//^D  i 


/ 

/  ®E  r- 


SOLID 

pg  -  Density 

C])  -  Compressional 
Wave  Speed 

Cg  -  Shear  Wave  Speed 


FIGURE  11.  A  COMPRESSIONAL  PLANE  WAVE,  k^,  IN  A  FLUID  HALF  SPACE 
IS  REFLECTED  FROM  A  FLUID-SOLID  INTERFACE  INTO  THE 
DIRECTION  kg  AND  REFRACTED  INTO  THE  SOLID  HALF  SPACE 
IN  THE  FORM  OF  A  COMPRESSIONAL  WAVE,  kg,  AND  A  SHEAR 
WAVE,  kg 
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The  incident  angle,  and  reflection  angle,  0g,  are  equal  while 
the  refraction  angles  are  found  from  a  Snell’s  Law  relation  to  be: 


and 


®D  = 

F 

o  •  -1  r^E  •  n  1 
Bg  =  sin  sin  6^J  . 


(76) 


(77) 


The  reflection  coefficient,  a,  in  equation  (74)  is  a  plane  wave 
approximation  and  hence  will  fit  nicely  into  the  empirical  model  as  a 
surface  distribution  coefficient  to  be  used  with  the  inclination  factor. 
This  will  permit  a  second  level  approximation  for  scattering  from  a 
finite  cylinder  providing  the  form  of  equation  (74)  can  be  modified  to 
correspond  to  the  parameters  used  for  expressing  the  structure  factor, 
equation  (40) .  The  only  parameter  i.i  equation  (74)  which  must  be  expressed 
in  wave-space  and  real-space  parameters  is  6^,  the  incidence  angle, 
since  0g  =  0^  and  0^  and  0g  may  be  derived  via  equations  (76)  and  (77). 

The  angle  of  incidence,  for  a  specific  differential  area  element,  can  be 
found  by  considering  the  inclination  factor  and  by  referring  to  Figure  4. 


cos  0. 
A 


(78) 


where  0^  is  the  angle  of  incidence  and  n  is  the  normal  to  the  differen¬ 
tial  area  element.  Equation  (78)  applies  equally  to  the  hemi-cylinder 
surface  as  well  as  the  end  disc.  The  form  of  a  is  different,  however, 
for  the  end  disc  versus  the  hemi-cylinder  when  expressed  in  wave  space 
and  real  space  parameters  because  of  the  difference  in  orientation  of 
the  surfaces  with  respect  to  a  given  ^  direction. 
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END  DISC  REFLECTION  COEFFICIENT 


Referring  to  Figure  5,  the  end  disc,  is  a  planar  surface. 
Thus,  the  unit  normals  for  all  differential  area  elements  point  in  the 
same  direction.  This  means  that  the  reflection  coefficient,  equation 
(7A) ,  is  independent  of  the  integration  variable  and  can  be  taken  outside 
the  integral  of  equation  (40),  i.e.; 


F«  (S)  =  a,  f  ^  dA,  (t)  .  (79) 

1  \i\  1 


The  form  of  for  the  end  plate  will  be  derived  for  the  back- 
scattering  case  only.  The  bistatic  case  has  been  investigated  qualita¬ 
tively  and  appears  to  be  different.  The  implication  of  this  difference 
appears  contrary  to  the  principle  of  reciprocity  but  is  still  under 
consideration. 


For  the  end  disc: 


cos  0, 


(80) 


Recalling  the  restriction,  described  earlier,  that  both  k^  and  k2 
are  coplanar  with  the  cylinder  axis  (hence,  they  define  the  k  -k  plane), 

u  A 

we  know  that: 


\t\^  =  k2  +  k2  .  (81) 

Therefore, 

sin  0.  =  [1  -  cos^  0.  ]”^ 

^1  n 
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=  [1  -  (7^)"]'' 

Ikl 

(82) 

Equation  (82)  can  be  used  with  equations  (76)  and  (77)  and  familiar 
trigonometric  double  angle  formulae  to  specify  all  the  parameters  in 
equation  (40)  for  the  reflection  coefficient.  The  results  are  listed 
below. 


.pAjii 


=  Ps^dH  -  , 


CpUl 


(83) 

(84) 


Z„  =  PqC„[1  -  (  ^  . 

h  =  c,  \t\ 


Cpk  , 


(85) 


2C„k 
E  X 


sin2(2e.  )  =  (-^)2  [1  -  (-^)2]  , 


C„k 

E  X  ^2- 


1  Cp  |k|  C.  |kl 


(86) 


coz2(2ep  )  =  [1  -  2(-^)212  . 

*^1  C.  |«| 


(87) 


Equations  (83)-(87)  may  be  substituted  into  equation  (74)  to  produce  the 
reflection  coefficient,  for  the  cylinder  end  disc.  The  second  level 
approximation  for  the  end  disc  structure  factor  is  found  using  equation 
(50): 
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(k)  =  (ap(2na2e 


2 

1^1  ^ 


]  . 


(88) 


HEHI-CYLINDER  REFLECTION  COEFFICIENT 

The  case  of  the  hemi- cylinder  reflection  coefficient  is  more 
difficult  because  ^2  depends  upon  the  integration  variable.  Recalling 
the  cylinder  orientation  and  the  parameters  shown  in  Figures  5  and  7,  it 
is  obvious  that  the  inclination  factor  varies  along  the  arc  length 
curvature  but  not  in  the  direction  of  the  cylinder  or  z-axis.  Hence,  we 
pick  a  representative  point,  P,  on  the  hemi-cylinder  surface  to  be  in 
the  X-Y  plane  as  shown  in  Figure  12. 

The  angle  between  the  scattering  vector,  and  the  unit  normal,  n, 
as  specified  in  equation  (78),  is  6^  as  shown  in  Figure  12.  However, 
the  wave-space  and  real-space  parameters  k^,  k^,  x,  and  z  can  be  expressed 
directly  in  terms  of  the  integration  angle,  (}»,  and  the  angle  6^  in 
Figure  12. 


=  cos  6.  . 
A 


But  from  spherical  trigonometry  we  have  the  relation: 


(89) 


cos  6^  =  cos  (}i  cos  .  (90) 

From  Figure  12,  we  see  that: 

cos  (])  =  -  (90a) 

a 

(Text  Continued  on  Page  48) 
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and 


Thus, 


cos  0,  = 


1  1^1  • 


1^1 


(90b) 


cos  0.  =  — - 

''  1^1  a 


(91) 


Equation  (91)  may  be  used  with  equations  (76)  and  (77)  and  the  double 
angle  formulae  to  provide  expressions  similar  to  equations  (83) -(87)  for 
the  cylinder  side  surface.  The  results  are: 


PpCpIkli 


(92) 


^2  ^  ^  4  \t\ei 


k  X  j 


(93) 


k  X 


®2  ®  ®  'f  iJja 


-h 


(94) 


2C  k  X  C  k  X 

sin2(20„  )  =  C-j5)2  [1-C^)2]  •  {l-(A21i.(^)21j  ,  (95) 

‘■2  '-F  |k|a  ‘•F  |k|a 


and 


k  X 

X 


cos2(26„  )  =  !l-2(A2[l-(-;f)^l)" 
®2  ‘=F  iJla 


(96) 


for  the  curved  side  surface  of  the  cylinder. 
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HEMI-CYLINDER  STRUCTURE  FACTOR 

The  structure  factor  for  the  hemi-cylinder  surface,  S2,may  be 
written  by  modifying  equation  (55)  to  include  the  reflection  coefficient, 
02 (x),  for  the  side  surface: 


F  it) 


k  L 


Ik^x 


2k  L  sin  (-|-)  a  «,Cx)  xe 
2 


0  (a2-x2)^ 


(97) 


where  a2(^)  is  given  by  equation  (74)  with  the  aid  of  equations  (92)-(96). 
Equation  (97)  must  be  evaluated  by  numerical  integration. 

COMPLEX  VALUES  OF  THE  REFLECTION  COEFFICIENTS 

The  reflection  coefficients  and  a2(x)  in  equations  (88)  and 
(97),  in  general,  are  complex.  For  incident  angles  0^,  in  equations  (76) 
and  (77),  which  are  less  than  the  critical  angle  for  compressional  and 
shear  waves,  i.e. , 


sin  [tt-^  sin  9.]  <  1  (98) 

^F,E  ^ 

the  reflection  coefficients  are  real.  For  6^^  greater  than  either  critical 
angle,  the  reflection  coefficients  are  complex;  i.e.,  real  plus  imaginary. 
In  order  to  evaluate  the  reflection  coefficients  in  these  complex  domains, 
the  refraction  angles,  6^  and  9j.,  were  expressed  in  terms  of  complex 
components : 


0.  =  Y.  +  i  p.  ,  j  =  B  or  E. 
J  J 


(99) 
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Thus,  the  Snell's  law  relations,  equations  (76)  and  (77),  take  the  form: 

C. 


sin  0.  =  sin  y-  =  sin  6. 

J  J  « 

(100) 

(for  Yj  <  n/2;  =  0) 

and 

C. 

sin  0.  =  cosh  p.  =-i  sin  0. 
j  Cg  A 

(101) 

(for  y.  =  n/2]  p.  >  0)  . 

J  o 

The  relations  to  be  used  in  the  reflection  coefficient 

for  0^  greater 

than  the  critical  angle  are: 

Z.  =  -iPjCjCsinh  , 

(102) 

sin^(20.)  =  -  sinh^(2pj  , 

J  J 

(103) 

and 

cos2(20.)  =  cosh2(2p.)  , 

J  J 

(104) 

where 

,  ,  C.  , 

P^  =  -jcosh  ^  [(^)  sin  0^]j 

(105) 

The  reflection  coefficient,  a,  in  equation  (74)  depends  upon  the 

densities  of  the  fluid  and  solid  (pj.  and  Pg) ,  the  compressional  wave 

speed  in  the  fluid  (Cj.) ,  and  the  compressional  and  shear  wave  speeds  (C^ 

and  Cg)  in  the  solid,  as  well  as  the  angle  of  incidence,  The  moduli 

and  phase  angles  of  the  reflection  coefficients  for  brass  and  aluminum 

planar  surfaces  in  water  are  shown  in  Figures  13  and  14,  respectively. 

The  modulus  and  phase  for  an  aluminum  surface  in  air  are  shown  in  Figure  15. 

(Text  Continued  on  Page  54) 
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(a)  Modulus 


Incident  Angle  (6^)  in  Degrees 


(b)  Phase  Angle 


Indicent  Angle  (0^)  in  Degrees 


FIGURE  14.  REFLECTION  COEFFICIENT  (a)  MODULUS,  AND  (b) 
PHASE  ANGLE  FOR  ALUMINUM  IN  WATER 
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The  moduli  for  brass  and  aluminum  in  water  are  similar  in  character. 
Each  has  a  small  cusp  at  the  critical  angle  for  compressional  waves  in 
the  solid,  followed  by  a  shallow  trough  which  rises  to  unity  at  the  cri¬ 
tical  angle  for  shear  waves  in  the  solid.  Phase  angle  plots  are  for  the 
absolute  value  of  the  phase  angle.  Hence,  the  large  cusps  in  the  phase 
angle  plots  are  artifacts  which  arise  from  the  folding  of  values  in  the 
range  180  degrees  -  360  degrees  back  into  the  range  180  degrees  -  0 
degrees  about  the  180  degree  line.  The  smaller  bumps  in  the  phase  angle 
plots  occur  at  the  compressional  critical  angle  while  the  nearly  360 
degree  phase  shifts  (large  artifact  cusps)  occur  at  the  critical  angle 
for  shear  waves. 

The  aluminum  plane  in  air.  Figure  15,  is  characterized  by  a  reflec¬ 
tion  coefficient  whose  modulus  is  unity  for  all  incidence  angles  and 
whose  phase  is  essentially  zero  everywhere.  This  could  be  anticipated 
since  the  ratio  of  the  densities  of  the  two  media  is  so  large.  Thus, 
the  case  of  aluminum  in  air  should  approximate,  very  closely,  a  rigid 
scatterer.  Indeed,  this  case  was  used  to  test  the  algorithm  for  com¬ 
puting  the  scattering  from  an  elastic  cylinder  via  the  reflection 
coefficients  as  outlined  previously.  The  results  are  presented  in  the 
next  section. 
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SECTION  VIII 

IMPEDANCE  CYLINDER  SCATTERING  CROSS  SECTION 


The  boundary  conditions  discussed  in  Section  VII  can  be  expected  to 
affect  the  dimensionless  scattering  cross-section,  Q(0 >  a  finite 
cylinder  providing  the  densities  and  wave  speeds  in  the  two  media  are 
sufficiently  different.  Variations  in  ft(0  due  to  the  scattering  ampli¬ 
tude  can  be  expected  to  be  moderate  when  no  dramatic  variation  is  observed 
in  the  modulus  of  the  reflection  coefficient  (See  Figures  13a  and  14a 
for  brass  and  aluminum  in  water).  However,  the  phase  variations  which 
occur  at  the  critical  angles  for  compressional  and  shear  waves  can  be 
expected  to  have  a  dramatic  impact  since  they  affect  the  manner  in  which 
two  or  more  scattered  waves  interfere.  This  hypothesis  suggests  that 
the  shear  critical  angle  will  be  the  dominant  parameter  because  of  the 
abrupt  2n  phase  shift.  This  is  consistent  with  Hickling's^  results 
which  indicate  that  the  spacing  between  peaks  in  the  experimental  fre¬ 
quency  response  for  elastic  spheres  is  most  strongly  correlated  with  the 
shear  wave  speed. 

The  scattering  cross-section  for  a  finite,  elastic  cylinder  with  a 
length  to  radius  ratio  of  6  was  calculated  using  the  end  disc  and  hemi- 
cylinder  structure  factors  given  by  equations  (88)  and  (97).  Evaluation 
of  the  integral  in  equation  (97)  was  accomplished  by  numerical  methods. 

Three  examples  were  investigated.  First,  Q(0  was  calculated  for 
an  aluminum  cylinder  in  air.  The  parameters  were  the  same  as  those 
specified  in  Figure  15.  The  purpose  of  this  example  was  to  test  the 
computer  algorithm.  In  principal,  the  results  should  be  the  same  as 
those  for  the  rigid  cylinder.  Comparison  of  the  scattering  cross- 
section  for  the  aluminum  cylinder  in  air.  Figure  16,  with  that  of  the 
rigid  cylinder.  Figure  8,  indicates  that  this  is,  indeed,  the  case. 

(Text  Continued  on  Page  57) 
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The  scattering  cross-section  for  the  aluminum  cylinder  in  water  was 
calculated  using  the  parameters  given  in  Figure  14.  A  three-dimensional 
plot  of  is  shown  in  Figure  17.  A  second  plot,  shown  in  Figure  18, 
has  been  rotated  180“  from  that  shown  in  the  previous  figure.  A  black 
and  white  contour  plot  showing  all  four  quadrants  is  given  in  Figure  19. 
Several  significant  features  are  clearly  visible  in  Figure  17.  First,  a 
single  broad  trough  extends  from  =  10  outward  at  approximately  30 
degrees  from  the  axis.  This  corresponds  to  the  critical  angle  for 
shear  waves  with  respect  to  the  side  surface  of  the  cylinder.  The 
breadth  of  the  trough  is  believed  to  be  due  to  the  variation  of  the 
shear  critical  angle  caused  by  the  curvature  of  the  side  surface. 

A  much  more  subtle  effect  is  observed  for  the  cylinder  end  disc.  A 

faint  ridgeline  can  be  seen  running  from  near  the  origin  diagonally 

across  the  maximum  which  is  parallel  to  the  t  axis.  This  ridgeline  is 

approximately  15  degrees  from  the  t  axis.  It  arises  from  the  small, 

z 

sharp  peak  in  the  reflectivity  curve,  Figure  14a,  which  corresponds  to 
the  critical  angle  for  compressional  waves  in  aluminum.  No  broad  minimum 
at  the  critical  angle  (30  degrees)  for  shear  waves  into  the  end  disc  is 
apparent.  Presumably  this  is  absent  for  two  reasons:  (1)  the  critical 
angle  is  sharp;  (2)  all  points  on  the  surface  of  the  end  disc  pass 
through  the  critical  angle  simultaneously  so  that  no  new  destructive  or 
constructive  interference  occurs  between  various  parts  of  the  end  disc. 
Since  the  critical  angles  for  the  end  disc  and  the  curved  side  surface 
do  not  overlap,  no  significant  cross  interference  between  the  contribu¬ 
tions  of  these  two  components  is  manifested.  That  is,  because  the 
amplitude  of  one  component  is  small  in  the  vicinity  of  the  shear  critical 
angle  of  the  other  component,  intensity  fluctuations  caused  by  cross 
interference  are  very  small. 

Another  subtle  feature  in  the  scattering  cross-section  for  aluminum 
in  water  is  the  shallow,  broad  minimum  in  the  ridgeline  along  the 

(Text  Continued  on  Page  61) 
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FIGURE  17.  SCATTERING  CROSS-SECTION  FOR  AN 
ALUMINUM  CYLINDER  IN  WATER 
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FIGURE  19.  SHADED  CONTOUR  PLOT  OF  THE  SCATTERING  CROSS' 
SECTION  FOR  AN  ALUMINUM  CYLINDER  IN  WATER 
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axis.  It  is  believed  that  this  minimum  is  coupled  to  the  flat  region  at 
30  degrees  arising  from  shear  wave  effects  and  that  it  will  show  oscilla¬ 
tory  behavior  with  c  complementing  flat  region  when  the  scattering 
cross-section  is  computed  to  larger  values  for  than  those  shown  in 
Figure  17.  One  final  point  can  be  made.  Since  the  approximation  utilized 
for  these  calculations  does  not  take  into  account  energy  which  penetrates 
the  surface  and  is  scattered  out  again,  the  region  in  Figure  17  from  30 
degrees  to  60  degrees,  measured  from  the  axis,  is  beyond  the  critical 
angles  for  both  the  end  disc  and  all  portions  of  the  curved  side  surface. 
Thus,  the  response  in  this  region  should  be  the  same  as  that  for  the 
rigid  cylinder.  Comparison  of  Figures  17  and  8  confirms  this.  Inspec¬ 
tion  of  Figures  18  and  19  reveals  no  additional  features  of  interest  for 
the  aluminum  cylinder  in  water.  However,  these  additional  respresenta- 
tions  do  help  to  clarify  points  already  described. 

The  scattering  cross-section  for  a  brass  cylinder  in  water  is 
presented  in  Figures  20,  21,  and  22.  The  parameters  used  in  the  evalua¬ 
tion  were  the  same  as  those  shown  in  Figure  13a.  In  contrast  to  aluminum, 
the  shear  critical  angle  for  brass  is  approximately  45  degrees.  Hence 
one  might  anticipate  that  there  would  be  strong  cross  interference 

effects.  The  spotty  minimum  which  runs  at  45°  to  the  ^  axis  in  Figure  20 

z 

is  believed  to  arise  from  this  cross  interference  between  the  end  disc 
and  hemi-cylinder. 

There  is  only  a  faint  hint  of  a  ridgeline  arising  from  the  peak  in 
the  reflectivity  curve  at  21  degrees  to  the  end  disc  near  the  origin  in 
Figure  20.  This  diminished  effp.ct  is  in  agreement  with  the  difference 
in  the  leflectivity  peaks  in  Figures  13a  and  14a.  The  beginning  of  a 
broad,  fiat  minimum  at  the  bottom  right  of  Figure  20  also  is  apparent. 

It  is  believed  that  this  is  the  second  oscillation  of  shear  effects 
mentioned  earlier  with  regard  to  Figure  17. 

Although  the  scattering  cross-sections  presented  in  this  section 
are  not  expected  to  describe  experimental  data  accurately,  they  do 
represent  an  improvement  over  the  rigid  cylinder  approximation. 

(Text  Continued  on  Page  65) 
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Additionally,  they  dramatize  the  significance  of  the  role  of  the  boundary 
conditions  in  determining  the  nature  of  the  scattering  cross-section  of 
any  target.  This  may  be  fortuitous  if  a  simple  set  of  target  indicators 
is  sought  which  emphasizes  universal  boundary  effects.  On  the  other 
hand,  if  geometry  indicators  are  sought,  it  may  represent  a  severe 
complexity. 
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SECTION  IX 

SUMMARY  AND  CONCLUSIONS 


Starting  from  the  wave  equation,  an  empirical  plane  wave  model  for 
acoustic  target  scattering  has  been  developed.  The  model  employs  plane 
waves  and  represents  the  total  scattered  field  as  a  sum  of  contributing 
effects  called  structure  factors.  The  principal  benefit  of  this  model 
is  that  it  provides  a  formalism  for  simultaneously  representing  both  the 
frequency  and  the  aspect  dependence  of  the  target  response  in  terms  of  a 
vector  parameter  in  wave-space.  The  relationships  between  the  wave-space 
and  real-space  parameters  have  been  described  so  that  features  which  are 
observed  in  the  target's  scattering  cross-section  in  wave-space  can  be 
related  to  real-space  parameters  such  as  target  geometry,  orientation 
and  boundary  conditions.  This  formalism  can  be  utilized  not  only  for 
further  refinements  of  the  empirical  model  but  also  as  a  means  for 
graphically  displaying  experimental  data.  In  this  way,  it  is  hoped  that 
interpretation  of  experimental  measurements  will  be  facilitated  by 
familiarity  with  similar  displays  of  model  results. 

The  model  was  applied  to  scattering  from  a  right  circular  cylinder 
of  finite  length.  This  cylinder  is  of  interest  because  it  is  a  mine-like 
target  and  because  it  has  a  non-analytic  surface;  that  is,  it  has  sharp 
edges.  It  is  these  edges  which  make  empirical  modeling  a  viable  approach 
to  the  acoustic  scattering  problem  because  sophisticated  theoretical 
models  break  down  when  mathematical  discontinuities  are  encountered.  A 
concise  mathematical  relation  was  derived  for  the  structure  factor  for 
the  right  circular  cylinder  with  rigid  boundary  conditions.  Elastic 
boundary  conditions  were  examined  and  incorporated  into  the  model. 

Because  of  the  complexity  of  the  mathematics  for  the  elastic  case,  a 
concise  mathematical  expression  was  not  possible.  Evaluation  of  the 
structure  factor  for  the  cylinder  with  elastic  boundary  conditions 
required  a  numerical  integration  by  computer.  One  effect  was  to 
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increase  the  computer  computation  time  for  evaluation  of  the  scattering 
cross-section  by  a  factor  of  32. 

Comparison  of  the  scattering  cross-section  for  the  cylinder  with 
different  boundary  conditions  indicated  dramatic  differences.  Although 
no  account  yet  has  been  taken  of  waves  internally  scattered  within  the 
target,  the  results  to  date  indicate  that  variations  in  boundary  con¬ 
ditions  alone  can  produce  markedly  different  scattering  results  for  the 
same  target  shape.  When  internal  scattering  is  included,  it  is  antici¬ 
pated  that  these  differences  will  increase  rather  than  dissipate. 
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